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[*"«- ■ A new method is developed for solving the conformally invariant integrals that arise 



in conformal field theories with a boundary. The presence of a boundary makes pre- 
■ vious techniques for theories without a boundary less suitable. The method makes 

essential use of an invertible integral transform, related to the radon transform, in- 
volving integration over planes parallel to the boundary. For successful application 
of this method several nontrivial hypergeometric function relations are also derived. 

> 

X 

1 Introduction 

At a critical point most statistical mechanical systems are not only scale invariant but 
are also conformally invariant |T], [|. This principle has profound implications for calcu- 
lations of the correlation functions, critical exponents and universal amplitudes of such 
systems 0. In two dimensions, where the conformal group is infinite dimensional, multi- 
point correlation functions are more strongly constrained then in dimension d > 2, where 
the conformal group is finite. However, consideration of d > 2 is also important, particu- 
larly in the statistical mechanical context when d = 3. In the case of general d conformal 
invariance still provides quite powerful constraints. For example, in the infinite geometry 
M. d the forms of the two and three point functions of scalar fields in a conformal field theory 
are determined exactly (up to normalisation) by the restrictions of conformal invariance. 
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Cardy has shown how to generalise the principle of conformal invariance to the case 
of the semi-infinite geometry Rl, so that surface critical phenomena can be probed using 
these techniques |I|, |]. In it is only appropriate to have conformal invariance under 
conformal transformation which leave the boundary fixed. In this case the restrictions 
on the form of correlations functions are not as strong. In particular the form of the 
two point function of a scalar field in IR+ is restricted by conformal invariance only up to 
some function of a single conformally invariant variable This function must be then 
be determined for the particular theory under consideration. 

In this paper we outline a powerful method, which makes essential use of conformal 
invariance, for calculating the two point functions of scalar, vector and tensor fields of 
conformal field theories in the semi-infinite space R^_. In particular we give a prescription 
for treating the conformally invariant integrals that arise in a diagrammatic expansion of 
the theory. Techniques for handling such integrals have been developed for the infinite 
space R d , and have proven to be very useful |5], ||. However these techniques do not extend 
to and so this alternative technique is developed. 

2 Conformal Invariance 

A transformation of coordinates conformal transformation if it leaves the 

line element unchanged up to a local scale factor Q(x). That is 

dx^dx 9 ^ = Q(x)~ 2 dx^dx^ . (2.1) 

For the discussion of two point functions of fields in a conformal field theory we need to 
consider the effect of conformal transformations on these fields. If a field 0(x) transforms 
under the conformal group as 

0{x) -> 0\x g ) = VL{x) ri O(x) , (2.2) 

for some rj, then 0(x) is said to be a quasi primary scalar field with scale dimension t]. 
A quasi primary vector field V^(x) with scale dimension r] is one which transforms like 

V(x) -> V^x 9 ) = Q(xyn, a (x)V a (x) , (2.3) 

where 7Z^ a (x) = Q(x)dxf l /dx a . The transformation for quasi primary tensor fields follows 
analogously. We will restrict our attention to quasi primary fields in this paper. 
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In the semi-infinite space we define coordinates x M = (y, x) where y measures the 
perpendicular distance from the boundary, and Xj are coordinates in the d—1 dimensional 
hyperplanes parallel to the boundary. The two point functions of scalar operators are 
restricted by translational and rotational invariance in planes parallel to the boundary to 
be 

(Oi(x)0 2 {x')) =G(y,y',\x-x'\) , (2.4) 

and scale invariance further restricts the form of G to depend on two independent scale 
invariant variables s 2 /y 2 and s 2 /y' 2 , where s 2 = (x — x') 2 . This situation should be 
contrasted with the case of infinite space where it is not possible to construct a vari- 
able from two points which is invariant under all of scale, translational and rotational 
transformations. 

For two points in conformal invariance provides further restrictions. Under con- 
formal transformations which leave the boundary fixed 

2 s 2 y_ , _ y' (2 ,\ 

S ~* Q(x)Q{x') ' y ~* Q{x) ' y ~* Q{x') ' 1 ' 

so that only one independent conformally invariant variable can be constructed from two 
points 

t=W OT " 2 = f = tt?' (2 ' 6) 

where s 2 = (x — x') 2 + (y + y') 2 is the square of the distance along the path between x 
and the image point of x' . 

As a consequence, the correlation function of two quasi primary scalar fields may be 
written as 

t0 1 (x)0 2 (x')) = , \ , \ f(0 , (2.7) 

for some arbitrary function /(£)[] 

As an example we consider free scalar field theory, where the field 4>(x) satisfies Dirich- 
let or Neumann boundary conditions at y = 0. Then by the method of images the Green's 
function is simply 

<0(*)0(*')> = G,(x, x>) = A(-L. ± ^) = j^iyrjm , (2-8) 



lr The ^ — > and £ — * oo limiting behaviour of this function is fixed by the Operator Product and 
Boundary Operator Expansions . 
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where 



A= {d _ 2)Sd ' V<t> = ¥-1, U(0 = C ± (I + 0-"" ■ (2-9) 



In the above expression the upper (lower) sign corresponds to Neumann (Dirichlet) bound- 
ary conditions and the factor Sd = 2n^ d /T(^d) is the area of a unit hypersphere in d 
dimensions. 

In henceforth referred to as /, the form of the two point functions of scalar, vector 
and tensor fields was worked out in detail for the O(N) sigma model in both the e and 
large N expansions. These calculation were significantly simplified by the use of a new 
technique to solve the conformally invariant integrals on that naturally arise. In the 
next section this technique is discussed in detail. 

3 Parallel Transform Method 

We consider integrals of the form 

4yz Ay'z 

where conformal invariance restricts the form of the integral to be a function of £ only. 
This follows because under conformal transformations which leave the boundary fixed, 
the integration measure transforms like d d x — > Q(x)~ d d d x and the factor l/(2z) d — > 
Q(x) d /(2z) d so the local scaling factor cancels. 

Given functions f\ and we may solve integrals of this type indirectly by first inte- 
grating /(£) over hyperplanes parallel to the boundary^ 

(y-y'Y 



/d d - 1 x/(0 = (%y') A /(p), p = A ^ L i A = |(d-1), (3.2) 

which defines the function f(p) to be 



fO>) = Wxjj Q duu^fiu + p). (3.3) 



2 This is related to the Radon transformation of f(x) which is defined as the integral of f(x) over 
all possible hyperplanes in R d . Here we consider integrals over the subset of hyperplanes in which 
are parallel to the boundary. 
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The crucial point is that this defines an integral transform / — > / which is invertible. 
Thus /(£) can be retrieved from f(p) via 



oo 

-A-l 



/lfl ^I dw ^ + «' (3 - 4) 

The integral in the above formula is actually singular for values of A that we consider 
here, but the inversion formula may still be defined by analytic continuation in A from 
Ke(A) < 0. To verify that the transformation (|3.3|) is compatible with the inversion 
formula ( |3.4j ) it is sufficient to make use of the following relation involving generalised 
functions 

du^p-uf^ul- 1 =B(fi,X)p^ +x - 1 ~ r(-A)r(A)5(p) as p -> -A . (3.5) 
For the case d = 3 when A = 1 we use 

^-~5'( P ) as A-l, (3.6) 



to reduce the inversion formula ( |3.4|) to the simple form 



f(o = --no- (3.7) 

7T 

Now that this parallel transform has been defined it is possible to derive an integral 
relation for the transformed functions by integrating /(£) in ( |3-1| ) with respect to x so 
that 

Kp) = f°°dz ^A(P)W') p = , ft = • (3.8) 

Jo 2z Ayz Ay'z 

In order to solve integrals of this type we first change variables z = e 2e , y = e 26 * 1 and 
yi = e 26» 2 go th a t equation ( |3.8[ ) becomes 

/( sinh 2 ^ - B 2 )) = J°Jl9 A (sinh 2 (# - 9 1 )) f 2 (sinh 2 (# - 2 j) . (3.9) 

Now by taking the Fourier transform 

f(k) = / dfle ifc 7(sinh 2 fl) , (3.10) 

J — oo 

then by the convolution theorem the transformed integral relation ( |3.8|) becomes 

}(k) = ) 1 (k)) 2 (k). (3.11) 
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Thus we may solve integrals of the general type given in by this double integral 

transform method provided that it is possible to make the transforms — > f%(p) 

fi(k) for both the functions fx and fi and that the subsequent inverse transforms of 

the resulting function f(k) can be made. Of course the form of the functions j\ and 
/ 2 are crucial in order for this procedure to be successfully undertaken. For the typical 
cases which arise in the diagrammatic expansion of a conformal field theory this method 
has proven to be very successful, although the intermediate steps often involve nontrivial 
manipulations of hypergeometric functions. In the next section several examples which 
are likely to occur in calculations in conformal field theory are given to illustrate the 
method, and provide a table of transforms for future reference. 



4 Illustration of the Method 

For application to the calculation of two point functions in a conformal field theory we 
may use this method to solve the integrals over products of propagators that occur in 
a diagrammatic expansion of the theory. Therefore, by considering, for example, the 



Green's function of the free scalar field given in fl2.8|) we wish to solve integrals of the 
following type 

o r I I 

dz d^r——- r^-^r , (4.1) 

J (2zY (s 2 ) a (s 2 ) Q (s /2 ) Q (s' 2 ) a 



with 



s 2 = (x - if + (y- zf , I 2 = (x - if + (y + zf , 

~s' 2 = (x' - if + (y' - z) 2 , f 2 = (x' - if + (y' + zf 



For conformal invariance, following ( [2.5D , we must also require 

a + a + a' + a' + (3 = d. (4.2) 

This integral may be readily cast into the general form (|3.1|) , for which we should then 
take 

^® = (2?/) a +«e*(i+£> ' = (2y') a,+a ' e a ' (i + ' ^ 

Later in this section we will consider the more general integrals that arise in the discussion 
of the large iV expansion of the O(N) sigma model, where the propagator for the auxiliary 
field A has a more complicated functional form. 
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To solve the integral fl4.1| ) using the method of section |3] we first take the sequence 

of transforms / — > / — > / as defined in ( |3.3| ) and Q3.10D for functions of the form 
above. For simplicity we take 

1 



f(£) = - 

The first transform / — > f follows from standard references 



(4.4) 



f(p) 



7T 



du u 



A-l 



1 



r(A) Jo (u + p) a (l + u + p) a 

A T(a + a — A) 



7T 



F^a + a — A, a ; a + a ; — 



(4.5) 



T(a + a) (l + p) Q +«- A " v ' l + p> 

The function F(a, b; c; z) is a hypergeometric function whose definition is given in ( |A.1|) 

For the subsequent transform / — > f we consider the cases a = 0, a = 0, a = a 
separately 



/i(0 

flip) 

?i(*0 
MO 

?n(*) 

/m(0 

/m(p) 



1 



7T 



Ar (a-A) 1 
r(a) p a " A 



(4.6) 



7T^ 



i d r(2a-2A)r(l-2a + 2A) 



r(a — Ah- 



+ 



r(a-\-±k) 



r(a)r(i + a-A) Lr(l-a + A + fA;) T(l - a + A - |fc) 



1 



(i + a 

X r(a-A) 1 



7T 



7T2 1 



r(a) (i + p)°- A 



(4.7) 



r(a)r(| + a - A) 



1 



r(a - A + £ife)r(a - A - I jfe) 



_ A r(2«-A) 



7T 



r(2a) (1 + p) 2q " a 



Ff2a — A, a ; 2a; 



1+P 



(4i 



1 ' r(a) r(± + ±A-J*)r(i + §A + ±*) 
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There is one other case, a particular combination of two functions of the type ( O ), which 
is of interest 



MO 



2£ + l 



r , , aT(2« - A - 1) 1 / 1 X , . 

fMp) = **> I F(2 tt -A-l, tt -l;2 tt -2;— ) (4.9) 



r(2a-l) (1 + p) 

/ IV (A;) = ;rt d 4 



d.a-iirdd - a) r(a - |A - I - j*)r(a - f A - ± + j*) 



r(a) r(jA - Jfc)r(|A + Jfc) 

The last two cases, /m and fiy, are important because the more general case where a 
differs from a by any integer follows in a straightforward manner from them. However, the 
derivation of those two results directly is nontrivial. The simplest way to verify them is 
by working backwards and taking the inverse transforms. A general procedure for taking 
the inverse transforms is discussed next. 



For application to conformal field theory where we have integrals of the form ( |3.1| ) 
then the transformed relation ( |3.11| ) suggests that we need to take the inverse transform 

of products of the functions fi(k) in I to IV. In all of these cases the dependence of f(k) on 
k is through combinations of Gamma functions. Consequently, by considering the poles 

of the Gamma function, the inverse Fourier transform / — > / of ( [3.11| ) can be performed 
by contour integration. We first consider the following combination of Gamma functions 
which is appropriate for verifying the transforms of fm and fiy above 

~ r(a-ifc)r(q + ifc) 

9a ^ k) ~ T{b-{k)V{b+\k) ■ (4 ' 1U) 

The poles of T(a — ~k) occur at ~k = a + n with residue ( — l) n /n\ (for n a non negative 

integer). Therefore, the inverse transform is obtained as a sum of the residues of g ab , 
resulting in a series that has hypergeometric form 

Msinh 2 #) = ^ fdke^ ke ~g afi (k) 

e - iaW F{2a,a-b+ 1; a + b; e~ m ) . (4.11) 



T(b-a)T(b + a) 

4r ( 2 «) 1 ^ „ , , l.o„ , o, n . 1 



F(2a,a + b- \;2a + 2b- 1; 



T(b ~a)T(b + a) {A cosh 2 6) 2a v ' 2 >""""" 'cosh 2 ; ' 

By choosing appropriate values for a, b, and noting that cosh 2 6 = 1 + p then the Fourier 
transformed functions / m and / IV follow directly from this result. To obtain the inverse 
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parallel transform we use 



1 f°° \ i 1 T(p + X) 1 . 



r(-A)y r " (i + p + O p r(p) (i + p+A 

with p = 2a + n, in the last line of Q4.11J ) so that 



42»-i7r A r(6 - a)r(6 + a) (l + £) 2a+A V 2 ' ' 1 + f ' 

(4.13) 

Now, with the appropriate choice of a, b, we can use this result to verify the parallel 
transforms fm and /rv in equations ( fO| ) and (^.9|). 



In order to solve the integrals of the type fl3.1|) we must find the inverse Fourier trans- 
form of products of the functions fi(k) in I to IV. These may can be simply obtained 
as hypergeometric series by contour integration in a similar way to above above calcula- 
tion. The procedure for finding the inverse parallel transform differs, though, because it is 
not always possible to make the simplifying manipulation of the hypergeometric function 
that is made in ( |4.11| ). This is because the hypergeometric series is often of higher order. 

However, a procedure for taking the inverse transform / — > / which bypasses this step is 
derived in the appendix. This procedure makes essential use of a special property of the 
hypergeometric series which arises on taking the inverse Fourier transform, that is due 

the symmetry g(k) = g(—k). After taking the inverse Fourier transform of products of 
the functions in I to IV, we obtain a hypergeometric series with one of the two following 
forms 

g(smh 2 9) = e- Aa W q+1 F q (2a,b 1 ,---b q ;c 1 ,---c g ;e-^), (4.14) 
/i(sinh 2 #) = e- 2a ^ q+1 F q (2a,b 1 ,---b q ;c 1 ,---c q] e~ 2 ^), (4.15) 

where the notation q +\F q refers to a generalised hypergeometric series which is defined 
in (|A.2[) . The crucial point is that the parameters hi and Cj in these functions are always 
related by q = 1 + 2a — b^ 

We now present the inverse transforms of six of the possible combinations of the func- 
tions in I to IV, which have been obtained using this method. These represent solutions 
to particular integrals of the type Q3.1]) . First we consider products of the functions f\ 
and /n- In these cases the inverse Fourier transform results in hypergeometric series of 
the form ( |4.15| ) and the inverse parallel transform can be found via the methods outlined 
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in the appendix. Thus, using ( |A.22|) , we obtain 



roc p 

= J dz Jd^v 



1 1 1 



J (2z) d I'af 

i jV(1 + a + a' — d)T(^d — a — a') , . , , 

^ - — r( i - — 1 * 1 +a+a ~ ft -« 

i d r(a + a — |c?)r(|d — a)r(|<i — a') 1 
+ 7F2 r(d-a-a')r(a)r(a') 

x F{\d - a, |d - a', 1 + |d - a - a'; -0 , (4.16) 
1 1 1 



/"OO /* 

j I)n (e) = y o d* ya*-^ 



id r(l + a + Q/-d)r(|d-a) 

r(|d)r(i + a'-id) ^^'"'2^ sJ» lLJ " 



/"OO P 

1u,u{i) = J dz jd d - l r 



1 1 1 

(2z) d (1 + (1 + £')«' 



= ^ + a + ^ ^ F(a, a'; 1 + a + a' - |d; -g) ■ (4.18) 
r(l + a + a' — |a) 

In order to bring these results to this form it is necessary to use several identities of the 
hypergeometric function which can be found in the standard references ||. 

If we take the limit a + a' — > d in these integrals, which corresponds to (3 — > in the 
original integral (|4.1| ) then the following relation 

' ^-S d 5 d (x-x') as /3->0, (4.19) 



can be used to show that 



jY(\d — a)T(^d — a') ,. . 

Xi,i + Xii,ii = tt" 12 2 >-5 d (x - aO , 4.20 

1 (a)L {a') 



in the limit a + a' —> d. This is the expected result when the range of the integral ( |4.1| ) , 
with a = at' = (3 = 0, is extended to the infinite space R d . In a similar way it is possible 
to show that if a + a' = d then X\^\ + Xn,i = 0, where Xn,i is defined by taking a <-> a 1 in 
Xyi. 

We now evaluate three more conformally invariant integrals involving combinations 
of the functions fm and fiy. In these cases the inverse Fourier transform results in a 
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hypergeometric series of the form ( |4.14| ). One obtains 

V{\d - a')F(a' - a)F(a + a' - A) i dAa >i d 1 



7r 2 d 4 a -5 d _ — (4.21) 



r(fd-a)r(a')r(a + !) [£(i + 0] 

'■ _|_ I 1 _|_ n, _ n,'- _ 

4^(1 + 

1 1 2f ' + 1 



x 3 F 2 (a, 1 + a — §d, |d — a'; a + |, 1 + a — a'; — ^ — -rj + a <-> a' , 



r(|d-a)r(|d-a / )r(a + a / -|d) irf 1 



r(a)r(a')r(rf-a-a') [£(1 + £)] 

x f(A - a, ~cf - a' ; d - a - a' ; -4£(1 + £) 

2f +1 2£' + 1 



^ 5d 777^-^TT7-I^ (4-22) 



2iv,iv(0 = fdz/d^ * * 



(2z) d £«(l + £)«£'«'(1 
r(|d - aQiy - a)r(a + a' - 1 - A) U 2£ + 1 

r(±d-a)rV)r(a-±) m + 0] a { ' } 

x 3-^2 (a, 1 + a — \d, \d — a'; a — ~, 1 + a — a'; ~ ^7^ + |j ) + 

To solve for Xni,m we require the transformed function fm(k) with the result for the 
inverse transform of the general case (|A.ll ) which is given in the appendix. For Xiv,iv 



we use fxy(k) with the inverse transform ( A..17 ). To obtain Xmrv i n the form ( [4.22 



we follow a similar procedure to the other two cases, but also use a relationship between 
hypergeometric functions with argument — z and hypergeometric functions with argument 
—1/z to simplify the expression. 

The solution to the integrals in (|4.16|) to ( |4.23| ) all have a pole at a = d/2 except 



for (|4.18| ). This pole arises due to the short distance logarithmic singularity for r ~ x in 



each of these integrals when a = d/2. 

We are now in a position to evaluate integrals of the type (|3.1j ) with products of more 
general functions than those discussed thus far. For example, if we consider the function 
g aj i given in fl4.13| ) which was derived from the definition of g ab in ( |4.1L1| ), then since 



~9 a ,b(k)~9 b ,c( k ) = ~9aM > ( 4 -24) 
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it follows directly that 

roo r I _ _ 

^ dzjd d -h—g a , b (0g b AO = 9a,c(0, a^c 

= (Ayy')^ d 5 d (x - x') , a = c. (4.25) 

This is a solution to an integral of the product of two hypergeometric functions with the 
special form ( |4.13|) . This relation is useful in the large N expansion of the O(N) sigma 
model with the Ordinary transition, where the Green's function of the auxiliary field A is 
a hypergeometric function of exactly this type JF), p| . 

We may generalise this further by considering the function 

* (k] = T(a-ik)T(a + {k)r(b-ik)T(b+ik) 

■ ' ~ r(c - \k)Y{c + ik)T(S - {k)T(S + \k)- 1 ' ] 

The methods of the appendix can then be used to obtain the inverse transforms of this 
function provided <5 = |Aor5 = | + |A. The inverse Fourier transform gives^ 

9abAsmh 9) ~ Y(c-a)Y(c + a)Y(6-a)Y(8 + a) e (4 ' 27) 
x 4 F 3 ( 2a, 6 + a, 1 + a — c, 1 + a — 5;1 + a — b, c + a, 5 + a; e -4 ' 6 *' 
+ a «-» 6 . 

Subsequently, using (|A.11| ) for the case 5 = ~ + |A we find 



„ ( ^ _J r(2a + A)r(&-a)r(6 + q) 1 ^ 

x 3 F 2 (a + |A, I + a - |A, c - b ; 1 + a - 6, a + c ; ~^j^y 
+ a <-> 6 , 

whereas when 8 — 5 A, using (|A.17) we obtain 



, £ x = 1 r(2a + A)r(6-a)r(6 + a) £+j , . 

</at,rfW - A^^Y(c-a)Y(c + a)Y(\\-a)Y(\\ + a) [£(1 + 

x 3^(1 + a + |A, 1 + a - |A,c - 6; 1 + a - 6,0 + 0;- ,^ ) 
+ a <->• 6 . 



3 This function, is related to the Meijer's G-function which is defined by the contour integration of 
combinations of Gamma functions with arguments of a particular form M . 
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Thus provided 5 is one of |A or | + §A then g a b,cs(0 can be obtained as 3F2 hypergeometric 
functions. The solutions to the integrals in (|4.21|) - (|4.23|) represent special cases of these 
functions. More generally, integrals of products of these types of 3 F 2 hypergeometric 
functions are possible. Since 

9(k) a b,cS9ce,bf(k) = 9{k)ae,6f , (4-30) 

then it follows that 

^ dz Jd d -h— d g abfi& (0g ceM (O = g ae , Sf (0 , (4.31) 

provided 5, f = ^A, |+|A. Similar integral relations can be derived by considering possible 
combinations of g a ^ with g a b, c s with particular choices of the parameters a, b, c, 5. Integrals 
such as these occur in a discussion of the large N expansion of the O(N) sigma model 
with the Special transition where the Green's function for the auxiliary field A contains 
hypergeometric functions of this type R, O, W2 . 



5 Integrals Involving Spin Factors 

We now turn our attention to conformally invariant integrals involving spin factors which 
occur in the discussion of two point functions of vector and tensor fields. For this we 

define fields X M , X^, with scale dimension zero, which transform like vectors at the point 
x under conformal transformations that leave the boundary fixed. 

X, = ,1*1 ^ > = JT^-JTtM ■ (5-1) 

£2(l+£)2 f 2 (l+0 2 

These are constructed to be unit vectors so that X^X^ = X^X^ = 1 . 

We will use the example of an integral with one spin factor in the integrand to illustrate 
the method. Such an integral would be appropriate for correlation functions involving a 
single vector field. We define 

POO f 1 _ ~ ~ 

I, = = I dz J d^r — X M A(0/ 2 (eO , (5-2) 

which has the functional form 1^ = J(£)X^ due to conformal invariance. To find /(£) we 
use the fact that X^ is a unit vector to obtain 

1(0 = fdz J d d -h ^L_(X • X)Mi)f 2 {?) . (5.3) 
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Now, since 



x . i= K+m+i)-&+i) t (5 . 4) 



4(«1+0((1+|) 



the methods of section |] can be used to solve for 7(£) in terms of hypergeometric functions. 
For example, if we take 

1 1 

then we may use the solution to the integral 2jri,rv m ( |4-22| ) to obtain 
7(0 = tt* 



d T(\ + \d-a)T(\d-a')T(\ + a + a' - \d) 1 



x 



r(| + a)T(a')T(\ +d-a-a') [f (i + 

F(d-2a,d- 2a'; \ + d - a - a'; -f) . (5.6) 

The case where a = ~(d — 1) can be worked out by an alternative method by noting that[| 

9, L.L r ^ = ^ d (* - , (5.7) 



W" 1 [|(1 + I)] i^ 1 ). 



so that from the definition of 7 M in Q5.2|) 



^((2^-i 7 N-^(2^[e(i+er ' (5i 



This result may be rewritten as a differential equation for 7(£) 

£(E(l+fl]* ( "''(fl)=*S,-— ^13, 



(5.9) 



which may be solved to give 



/(0 = , 1 ; — — F(l,d-2a';l + y-a'; -£) . (5.10) 

vs; d - 2a' + £)] a "l V 2 y 

The constant of integration is taken to be zero, because otherwise the presence of such a 
term would violate ( |5.8| ) by producing an extra delta function contribution to the RHS. 
This solution is in agreement with ( |5.6| ). A similar procedure can be used for integrals 
involving more spin factors, which would be appropriate for correlation function involving 
the energy momentum tensor or two vector fields for example. Such integrals are evaluated 
in 7 by a slightly different method. 

4 For this we recall that -d 2 s 2 ~ d = (d - 2)S d 8 d {s) . 
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6 Large N Expansion for the 0(N) Model 



In this section, by way of conclusion, we demonstrate the use of the of the parallel trans- 
form method to to calculate two point functions in the expansion of the O(N) non-linear 
sigma model for the case of semi-infinite geometry. As usual, the nonlinear constraint on 
the fields <p a (x)\ <\> 2 = N can be removed by introducing an auxiliary field X(x) in the 
Lagrangian via an interaction term C\ = ^X(p 2 . To analyse the two point functions of the 
fields 4> a and A we first define 

{<p a (x)(j)p(x')) = G^(x, x')<5 a /3 , (X(x)X(x')) = G x {x,x') . (6.11) 

Then, to zeroth order in the i expansion, these Green's functions satisfy the following 



relations [13 



( - V 2 + (A(x)))G^(x, x') = 5 d (x - x') , (6.12) 
d d rGl{x,r)G x {r,x') = -^S d (x - x') . (6.13) 

Both of these relations may be solved by making use of conformal invariance and using 
the parallel transform method discussed in section || For this we write 

G * (z ' x " ] = ' Gx(x - x ' ] = jdr^ m) ' (fU4) 



Since 2?7^ + rj\ = d due to conformal invariance of the integral in Q6.13 ), then the zeroth 
order result rj^, — \d — 1 implies that r] X = 2 to this order. Now, with the scaling relation 
(X(x)) = A\/4y 2 , it is possible to obtain as a solution to a differential equation. 
Alternatively we can recast ( |6.12j ) into an integral equation so that the method of parallel 
transforms can be used to obtain a solution. Writing 

d d r H{x, r)G^(r, x') = 5 d (x - x') , (6.15) 

requires that 

H(x,x') = (- V 2 + ^)5 d {x - x') . (6.16) 
The integral of H (x, x') over planes parallel to the boundary may be written as 

d d -^H(x,x') = —L-rhfat/) , (6.17) 
[Ayy') 2 
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defining h to be 



h(e™,e™') = (-^- 2+ l + A x )5(6-6>). 
The subsequent Fourier transform of h(e 2e ,e 28 ') gives the simple expression 

h(k) = k 2 + 1 + A x . 



(6.18) 



(6.19) 



We may now solve for G$ by first integrating the integral equation ( |6.15 ) over planes 
parallel to the boundary and then taking the Fourier transform as defined in ( 3.10| ). The 
resulting equation is 

h(k)Uk) = 1 , (6.20) 

where f^k) is the transform of the function /,/,(£) defined in ( |6.14| ). Consequently the 
desired result is 

1 1 



/*(*) 



h{k) 



k 2 + 1 + Ay 



(6.21) 



If we are now express fAk) as 



h(k) 



i r(// + |fc)r(//-|fc) 
wr(i + ii + ik)r(i + n-{k) 



l + A) 
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(6.22) 



then we may use the result ( I.13f) to obtain the inverse transform directly: 



i r(2/i + A) 



4 i+2 M7r A r (! + 2/i ) 



TT 1__F(2 / , + A, | + 2 M ; 1 + 4^; ^ 



(6.23) 



This general form for /</,(£) gives the correct large N Green's function G<p(x, x') appropriate 
for both the Ordinary and Special transitions in the statistical mechanical context where 
we should take \i — (d — 3)/4 and \i = (d— 5)/4 respectively |T^| . Solutions for G\(x, x') 
can now be obtained in a similar way via the integral equation Q6.13 ). Results for G\(x, x') 
for both the Ordinary and Special transitions were calculated with the parallel transform 
method in / and also in [11], [L2| by a different method. It would be interesting to see 
if the next order in the expansion can be obtained using the methods discussed in this 
paper; this is the subject of future research. 
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Appendix: Hypergeometric Function Relations 

In this appendix we derive some essential hypergeometric function relations that are need 
in section |j. We start with the definition of the hypergeometric function 

n^;c;^E^, (A.i) 

where (a) n = T(a + n)/T(a) is the Pochhammer symbol. There is a natural generalisation 
of this definition, which is called a generalised hypergeometric series 

t-t / \ Z ' ' ' ( a p)n / A o\ 

p F q (a u • • • a p ; ci, • • • c q ; z) = }^ r~? \ • (A.2) 

n=0 n - \ C l)n ' ' ' \ c qjn 

For application to section [| we need to consider the inverse parallel transform of functions 
of the following hypergeometric form 

£(sinh 2 #) = e- Aa ^ q+1 F q (2a,b 1 ,---b q ;c 1 ,---c q ,e-^), (A.3) 
h(smh 2 6) = e- 2a ^ q+1 F q (2a,b 1 ,---b q ;c 1 ,---c q ;e- 2 ^), (A.4) 

where Cj = 1 + 2a — hi. To take the inverse transform we express these results as the sum 
indicated by (|A.2|) and then observe that, with p = sinh 2 9, 



e -2p\e\ 



p + ^+~py 2p = ^TojA^ 5 + p; 1 + ^ > ( A - 5 ) 



where p = 2a + 2n for g and p = a + n for h. We can now obtain the inverse transform 
by using 

i ' = W' , (A ,) 



r(-A)A (i+/> + 0» r( P ) (i + 0" +A 

with the result 

r(/ ' } ^pp^ 1 - ■ 1 ■ F(p, k+P;l + 2p: 1 
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= m +^L»» F & + A + 1)1 1 + w " A); 1 + p • ~ (A ' 8) 
= ^rb^ F ^ p + 1(1 + p ~ A); 1 +p : • (A ' 9) 

The techniques for finding the inverse transform of g and h are related, except that for h 
it is necessary to use (|A.7| ) while for g either of the equivalent results ( |A.8| ), (|A.9|) may be 



used. Both of these equivalent results are helpful, because two different expressions for 
g(£) can be derived from them and a nice simplification of these expressions occurs for 
different values of the parameters b{. For the purpose of this discussion we will focus on 
the inversion of g for which we use ( |A.9|) with p = 2a + 2n to obtain 



1 r(2a + A) 1 ~ 1 (2a) n (fe 1 ) n -- ■ (b q ) n (2a + A) 2n 

_ " 4^ r(2o) [e(l + 0] a+ ^ A n=o^ (ci)n " " " (c q ) n 4«(2a) 2n 1 ' j 

X [^T^jp F ( Q + |A + n, i + a - |A + n ; 1 + 2a + 2n ; - — 
By expanding the hypergeometric function, g(£) can be rewritten as 

o(f) = 1 F ( 2a + A ) 1 V C (~ l)N (A 11) 

m) 4*^ r(2a) [^i+or^jfci, [4£(i+or 1 j 

where the coefficient Gat is given by the finite sum 

G " (-1)" (2a) n (b 1 ) n ---(b q ) n {2a + X) 2n {a+\X + n) N _ n {\ + a-\X + n) N _ n 
N ^ n\{N-n)\ (ci) n --.(c,) n 4«(2a) 2n (l + 2a + 2n) iV -n 

(A.12) 

This can be simplified further by using the following identities for the Pochhammer symbol 



l n (b)n(| + \p)n , (P + n) N -n = ^ , (p)-» = 



so that becomes 

GW = (Q + l2 Nul+ + 2a)~ ' A)jV ^ +3 (2a, 1 + a, b u - ■ ■ , b 9 , \ + a + ±A, -N; (A. 14) 

a, ci, • • • , c q , I + a — |A, 1 + 2a + N; 1 

The finite q+ ±F q+ 3 finite hypergeometric series, with argument 1, has a special form be- 
cause Ci = 1 + 2a — 6j. As a consequence, for the particular case g = 1, the resulting 5 F 4 
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can be summed exactly by a special limit of Dougall's theorem which states that |T5 



(1 + 2a) N (l + 2a - b - c) N 

5 F 4 2a, 1 + a, 6, c, -iV; a, l+2a-6, l + 2a-c, l + 2a + iV; 1 = ) ,N \ . 

(1 + 2a — o)jv(1 + 2a — c)n 

(A.15) 

So for q = 1 and applying Dougall's theorem we find 

1 T(2a + A) 1 „ / 1 , 1 1 , , , 1 

yv ^ ; 4 2a 7r A r(2a) [^+l)] a+ i A V 4^(1 + 

(A.16) 

This result and the inverse fourier transform ( 4.11|) are sufficient for verifying the integral 
transforms in ( |4.8| ) and (|4.9|) . 

Although a generalisation of Dougall's theorem is not known for arbitrary q, the 
coefficient Gn can usually be simplified to give a finite 5F4 series through cancellation of 
the parameters. However, for this cancellation to occur it may be necessary in some cases 
to use an alternative formula for g(£) which may be derived from (|A.8|) following a similar 
procedure. The result is 

_ 1 r(2a + A) £+J » (-1)" 

r(2a) [e (i+e# +o+ ^^a [4ea+or ' ( j 

where 

Giv = ^ + Q ^j"^ ° " g+ 4^ +3 (2a, 1 + a, 6 X , • • • , b q , a + |A, -N; (A. 18) 

a, ci, • • • , c q , 1 + a — |A, 1 + 2a + N ; l) . 



The calculation of a(£) proceeds in a similar way. Using ( |A.7 ) one gets 



1 I> + A) 1 » (-1)" 



where 



H N = {a ^+2^ (2a, 1 + a, 6 1; ■ • ■ , 6 g , -N; (A.20) 

a, ci, • • • , c q , 1 + 2a + N; lj . 
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Again we recall that Cj = 1 + 2a — b^. For the case q = 1 which is relevant for section ^ 
we use a theorem similar to Dougall's, 



4 F 3 (2a, 1 + a, 6, -JV; a, 1 + 2a - b, 1 + 2a + JV; 1 



) 



(l + 2a)jv(| + a-6) JV 
(i + a)^! + 2a - 6)jv ' 



(A.21) 



to obtain 



HO 



1 T(a + A) 1 



F(a + A, ± + a-o;l + 2a -&;--). (A.22) 



4 a 7r A T(a) £ a + 
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